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ABSTRACT 

This paper presents late-time near-infrared and Spitzer mid-infrared pho- 
tometric and spectroscopic observations of warm dust in the Type Iln SN 
2005ip in NGC 2906. The spectra show evidence for two dust components 
with different temperatures. Spanning the peak of the thermal emission, these 
observations provide strong constraints on the dust mass, temperature, and 
luminosity, which serve as critical diagnostics for disentangling the origin and 
heating mechanism of each component. The results suggest the warmer dust has 
a mass of ~ 5 x 10~ 4 M© and originates from newly formed dust in the ejecta, 
continuously heated by the circumstellar interaction. By contrast, the cooler 
component likely originates from a circumstellar shock echo that forms from the 
heating of a large, pre-existing dust shell ~ 0.01 — 0.05 M© by the late-time 
circumstellar interaction. The progenitor wind velocity derived from the blue 
edge of the He 1 1.083 pm P Cygni profile indicates a progenitor eruption likely 
formed this dust shell ~100 years prior to the supernova explosion, which is 
consistent with a Luminous Blue Variable (LBV) progenitor star. 
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1. Introduction 

SN 2005ip was discovered in NGC 2906 (d « 30 Mpc) on UT 2005 November 5 (Boles 
et al. 2005). Early optical spectra suggested the discovery occurred a few weeks following the 
explosion (Modjaz et al. 2005). The supernova is a Type Iln event given the development 
of narrow Ha emission lines (Smith et al. 2009b). Fox et al. (2009) discovered a late-time 
near-infrared excess that has persisted at least two and a half years post-discovery. 

Late-time infrared emission typically indicates the presence of warm dust. The origin 
and heating mechanism of the dust, however, is not always well constrained. The dust may be 
newly formed or may have existed at the time of the supernova. If newly formed, the dust may 
condense from refractory elements in the expanding supernova ejecta or in the cool, dense 
shell of post-shocked circumstellar gas lying in between the forward and reverse shocks (like 
the post-shocked environment trailing stellar wind collisions in WR binary systems (Usov 
1991)). In both cases, several heating mechanisms are possible, including radioactivity, 
optical emission from circumstellar interaction, and collisional heating by hot gas in the 
reverse shock. 

Alternatively, pre-existing dust may be eollisionally heated by hot, shocked gas or ra- 
diatively heated by either the peak supernova luminosity or late-time optical emission from 
circumstellar interaction. In the latter case, an ‘IR echo’ is evident due to light travel time 
effects. Multiple scenarios can contribute to the late-time infrared flux, as in the cases of 
SNe 2004et (Kotak et al. 2009, Sugerman et al. in prep) and 2006jc (Mattila et al. 2008). 

Type Iln supernovae, named for their "‘narrow” emission lines (Schlegel 1990), are more 
often observed to exhibit late time infrared emission associated with warm dust than any 
other supernova subclass (e.g., Pastorello et al. 2002; Gerardy et al. 2002; Pozzo et al. 2004; 
Fox et al. 2009; Smith et al. 2010; Miller et al. 2010a, b). Disentangling the origin and heating 
mechanism of the dust can yield important diagnostics concerning the circumstellar medium, 
supernova progenitor, and explosion dynamics. For example, if the dust is shock heated, the 
observed dust temperature yields the gas density (Dwek 1987; Dwek et al. 2008), which can 
be used to trace the progenitor’s mass loss history (e.g., Smith et al. 2009b). If the dust is 
newly formed, the observed dust mass can be compared to models that predict supernovae 
as primary sources of dust at high redshifts (Todini & Ferrara 2001; Nozawa et al. 2003, 
2008). 

The origin of the observed dust is not always obvious. For SN 2005ip, optical spectra 
show a progressive attenuation of the red wing of both the broad and intermediate lines, 
directly confirming the formation of new dust in both the ejecta and post-shocked cool, 
dense shell (Smith et al. 2009b). Such direct evidence, however, is rare. No more than a 



handful of supernovae (Lucy et al. 1991; Meikle et al. 1993; Elmhamdi et al. 2004, 2003; 
Sugerman et al. 2006; Kotak 2008; Pozzo et al. 2004; Smith et al. 2009b, 2008b) show direct 
evidence of dust formation in the ejecta, and aside from SN 2005ip, only SNe 1998S (Pozzo 
et al. 2004) and 2006jc (Smith et al. 2008b) show direct evidence of dust formation in the 
cool, dense shell (although this scenario is invoked to explain the asymmetries observed in 
the spectra of SNe 2004et, (Kotak et al. 2009) and 2007od (Andrews et al. 2010)). The 
observed dust yields, however, all tend to be 2-3 orders of magnitude lower than required to 
account for the large amounts of dust observed at high redshifts. 

Mid-infrared (3 < A < 15 pm) observations span the peak of the thermal spectral energy 
distribution from dust with temperatures ranging 100 < T d < 1000 K, providing strong 
constraints on the dust mass, temperature, and, thereby, the luminosity. These quantities 
serve as useful diagnostics for disentangling the origin and heating mechanism of warm dust. 
Late-time mid-infrared observations of supernovae, particularly the Type Iln subclass, are 
rare. 


In this paper, we present the first late-time (day 936 post-discovery) Spitzer / IRS spectra 
of a Type Iln supernova, as well as coincident Spitzer / IRAC photometry. We also present 
0.9-2. 5 pm, R = 3000 spectra obtained with APO 3.5-m/TripleSpec. Section 2 presents the 
observations and data reduction techniques. The combined spectra show evidence for two 
independent dust components: a hot, near-infrared (HNI) and warm, mid-infrared (WMI) 
component. For each component, we derive the dust composition, mass, temperature, and 
luminosity. Section 3 explores the origin and heating mechanism of these components to 
determine the degree to which SN 2005ip forms new dust. The HNI dust mass originates 
predominantly from newly formed dust in the ejecta, while the WMI component likely orig- 
inates from an ‘IR echo’ in a pre-existing dust shell. We use these results to explore the 
progenitor system and its evolution. Section 5 presents a summary of the findings and a 
discussion of future work. 


2. Observations 
2.1. Spitzer 

As part of PID 50256, the Spitzer Infrared Spectrograph (IRS) (Houck et al. 2004) 
obtained one mid-infrared spectra on June, 3 2008 (936 days post-discovery) with the Short- 
Low module (SL. R.~ 60-T20, 5.2-14 pm). The Spitzer Infrared Array Camera (IRAC) (Fazio 
et al. 2004) followed the IRS observations with images of the supernova and host galaxy in 
all four bands on June 10, 2008. Table 1 lists the observational details. 



2.1.1. IRAC 


Spitzer collected 5-minute IRAC exposures consisting of ten 30-second integrations. 
Pipeline-reduced and calibrated (BCD) images were taken from the Spitzer archive, and 
combined into single frames with enhanced pixel resolution of 0'. / 75 pix -1 using the MOPEX 
software package provided by the Spitzer Science Center. Figure 1 shows a post-BCD 3.5 
pm IRAC image with the IRS map overlay. Since the supernova lies near the center of 
the NGC 2906, the rapiclly-varying background of the host galaxy complicates aperture 
photometry. Instead, a number of unsaturated, linear, and isolated (e.g. no other sources 
in the wings) stars were used to build an empirical PSF for each channel, each of which 
was used to measure the brightness of the supernova. In each channel, the residuals from 
subtracting the best-fit PSF were small compared with the predicted uncertainty that the 
photometry task allstar provides, which factors in Poisson noise along with flat-field and 
profile-fitting errors as well as read noise. PSF-fit measurements of field stars in the frame 
were consistent with aperture photometry to within 5% in all channels. 


2.1.2. IRS 

The Spitzer / IRS mapped the position of SN 2005ip, with 12 cycles of 5 explosures, 
stepped 2" 7 perpendicular to the 3'.'7 wide slit. The target spectrum was extracted and cali- 
brated from the central pointing observation using the SMART data analysis software (Higdon 
et al. 2004), with bad pixels identified by IRSCLEAN. Subtracting off-order observations re- 
moved sky and zodiacal background. All data collection events for a given order were then 
median combined and the spectra extracted using the advanced optimal extraction routine 
AdOpt (Lebouteiller et al. 2009). The AdOpt tool provides the powerful ability to simulta- 
neously fit multiple sources using a super-sampled PSF plus a complex background at each 
row. This feature disentangles the supernova emission from that of the galaxy nucleus as 
well as removes background emission of the galaxy’s spiral arms (Figure 2). 

Figure 3 shows the resulting IRS spectrum redshifted to account for the radial velocity 
of the galaxy (2140 km s^ 1 -}- Vlsr ), along with the IRAC photometry. Continuum emission 
peaking at around 3-4 pm tends to dominate the spectrum. Few, if any, emission lines are 
apparent. 
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2.2. TripleSpec 

TripleSpec, an 0.9-2. 5 pm, R = 3000 spectrograph operating at APO (Wilson et al. 
2004; Herter et al. 2008) obtained a spectrum on day 862 post-discovery. Forty minutes 
of on-source integration consisted of 8 independent 5 minute exposures nodding between 2 
different slit positions. We extract the spectra with a modified version of the IDL-based 
SpexTool (Cushing et al. 2004). The underlying galactic arm and sky are approximated in 
SpexTool by a polynomial fit and subtracted from the supernova. TripleSpec observations 
on day 1243 post-discovery suggest little evolution occurred between the two epochs. We 
therefore create a single spectrum from the near- and mid-infrared spectra (see Figure 3). 


2.3. Dust Composition, Temperature, and Mass 


Assuming only thermal emission, the combined near- and mid-infrared spectra provide a 
strong constraint on the dust mass, temperature, and thereby, the luminosity. The luminosity 
of a single spherical dust particle of radius, a , and temperature, T d , is given as 

L d { A) = AKa\ixB u (T d )Q u {a)), (1) 


where B v (T d ) is the Planck blackbody function and Q„(a ) is the emission efficiency. Fox 
et al. (2009) and Smith et al. (2009b) provide two pieces of evidence for optically thin dust: 


1) T 


^NIR 


0.5 and 2) the relatively high transmission of X-rays responsible for 


^NIR+ioPT 

ionizing the unshocked circumstellar medium. For optically thin dust with mass, M d) at a 
distance, d, from the observer, thermally emitting at a single equilibrium temperature, the 
total flux can be written as: 

M d B v (T d )n v {a) 


F u 


d 2 


where K u {a). the dust mass absorption coefficient, is: 

3 




(2) 


( 3 ) 


for a dust bulk (volume) density p. 


Given simple dust populations of a single size composed entirely of either silicate or 
graphite, Figure 4 plots the dust absorption coefficient and emission efficiency for several 
grain sizes of each composition, which is derived from Mie theory. For the observed spectrum 
in Figure 3. Figure 5 shows the best fit of equation 2 with graphite and silicate models using 
IDL’s MPFIT function. The lack of an emission feature at ~9 pm immediately rules out any 
silicate grain contribution. We therefore use only graphite models throughout the rest of 
this paper. 
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We fit a multi-component model to the combined spectrum, where each component as- 
sumes a single dust mass, M d , at a single temperature, T d , composed of graphite of a single 
grain size. The dust mass, temperature, and grain size of each component are all free param- 
eters. Figure 5 shows the optimal fit, which consists of two dominant components: a “hot” 
(~800 K), near- infrared (HNI) and “warm” (~400 K), mid-infrared (WMI) component. 

Adding additional components to the fit tend to yield relatively small dust masses that 
do not improve x 2 by more than a couple percent. This result has two implications. First, 
two unique dust components exist, as opposed to a single dust component with a continuous 
temperature distribution. Second, a hotter (> 1000 K) third component, if it exists, does 
not significantly contaminate the near-infrared fits at these late times. (At early times, a hot 
photospheric component is expected and is, in fact, observed to dominate the early J-band 
data in Fox et al. (2009).) We therefore only consider two-component fits throughout the 
rest of this article. 

Table 2 lists the dust masses, temperatures, and luminosities associated with various 
grain size combinations, but overall, grain size has little consequence on either these parame- 
ter values or the resulting goodness-of-fit (x 2 ). The only size dependent variable in equation 
2 is k, but Figure 4 shows the dust opacity coefficient for graphite is independent of grain 
radius at infrared wavelengths (>1 gm) for grain sizes <1 gm, which is typical for most 
grains. Other methods (discussed in Sections 3.2 and 3.3), however, can constrain the grain 
size. 

Figure 6 plots the luminosities of both the HNI and WMI components on day 936 for 0.1 
/.mi grains, along with the time-series evolution of SN 2005ip at near-infrared (Fox et al. 2009) 
and visible luminosities (Smith et al. 2009b). Smith et al. (2009b) show that as the photo- 
spheric component drops off over the first ~100 days, an optical luminosity plateau, L p ] ateau , 
appears and extends throughout the extend of the observations. This plateau arises from 
radiation generated by continuous shock interaction with the dense circumstellar medium, 
as opposed to an optical light echo powered by the peak supernova luminosity. 

Lacking mid-infrared observations, the data obtained prior to day 936 cannot distin- 
guish between the multiple dust components. The HNI component in Figure 5, however, 
dominates the near-infrared observations (i.e., little contribution from the WMI component 
at near-infrared wavelengths), which plateau throughout the extent of the observations (see 
Figure 6). In fact, additional epochs of TripleSpec spectra obtained on days 862 and 895 
show little evidence for spectral evolution at later times. Given these observations, we as- 
sume throughout this paper that all near- infrared results predominantly represent the HNI 
component evolution. For both the HNI and WMI components, the thermal emission arises 
from warm dust. Disentangling the composition, origin, and heating mechanism of these 
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dust components, however, requires a detailed analysis of potential heating mechanisms. 


3. Analysis: Dust Origin and Heating Mechanism 

Section 1 summarized likely origins and heating mechanisms for late-time infrared dust 
emission, distinguishing scenarios involving newly formed versus pre-existing dust. Smith 
et al. (2009b) provide spectroscopic evidence for new dust condensation in the ejecta from 
days 60-170 and in the cool, dense shell of post-shocked gas at days >413. This dust forma- 
tion timeline closely corresponds to the evolution of the near-infrared light curve in Figure 
6. The relative contribution of this newly formed dust to the HNI and WMI components, 
however, is not immediately clear, but can be addressed by the Spitzer mid-infrared pho- 
tometry and spectroscopy presented here in the context of both pre-existing dust scenarios: 
shock heating and an infrared echo. 


3.1. The Blackbody and Shock Radii 

For both shock heating and an infrared echo, the blackbody and shock radii respectively 
serve as useful reference points. The blackbody radius, given as r bb = ( ) 5 where a is 
Stefan-Boltzmann’s constant, defines the minimum shell size of an observed dust component. 
In the case of SN 2005ip, blackbody fits ( Q = 1) of the combined spectrum on day 936 in 
Figure 5 yields blackbody radii of ry,(WMI) «4.8 x 10 16 cm (0.048 ly) and r^HNI) «7.7 x 
10 15 cm (0.0077 ly). 

The shock radius, given as r s = v s t for a constant velocity, v s , defines the maximum 
radius that the forward shock can travel in a time, t. The shock velocity can be derived from 
the optical emission line widths. For SN 2005ip, Smith et al. (2009b) observe broad emission 
lines corresponding to radial velocities of 15,000 km s 1 (0.05c) through 900 days post- 
discovery, yielding a maximum shock radius on day 936 of r s \ ~ 1.25 x 10 17 cm (0.125 ly). At 
the same time, the intermediate width lines correspond to slower shock velocities of ~1000 
km s^ 1 (0.003c) through ~900 days post-discovery. Chugai & Danziger (1994) propose that 
two unique shock velocities can coexist if the progenitor’s wind is clumpy or asymmetric, as 
opposed to homogeneous and spherical. A relatively rarefied wind allows uninhibited shocks 
to maintain the observed high velocities, while a much slower shock propagates through the 
denser regions. 

Since a shock will typically destroy any pre-existing dust, understanding the shock 
propagation is essential for modeling the evolution of any late-time infrared emission. If, 
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for example, the denser regions are concentrated in a homogeneous equatorial disk, only 
the slower shocks need be considered as they propagate contiguously throughout the disk. 
Alternatively, if the denser regions are concentrated in clumps, the situation is more compli- 
cated. The fastest shocks will reach the clumps first, at which point the shock velocity drops. 
The average shock velocity is a function of the clump filling factor and distribution. Since 
the distribution of the dense regions is not well known, we consider a second shock radius 
for which we adopt a more modest shock velocity of v s ~5000 km s~ l , which is consistent 
with the intermediate component of many Type Iln supernovae, (e.g. Salamanca et al. 2002; 
Smith et al. 2008a, 2010; Steele et al. 2008). This velocity yields an upper limit on a second 
shock radius on day 936 is r s2 ss 4 x 10 16 cm (0.04 ly). 


Both the HNI and WMI shells are optically thin, thereby confirming the assumptions 
presented in Section 2.3. The optical depth of each shell can be written as, 


M d 

- 47rr 2 K -g 


(4) 


where K avg = 435 cm 2 g~ x is the absorption coefficient for graphite averaged over 1 — 15 //in. 
For 0.1 pm grains, Table 2 shows M d (HNI ) ~ 5.2 x 10“ 4 M 0 and M d {WMI) « 4.3 x 
1CT 2 M 0 . Given the minimum radii, r > r^HNI ) «7.7 x 10 l£> cm (0.0077 ly) and r > 
rbb(WMI ) :=s4.8 x 10 16 cm (0.048 ly), equation 4 yields r < 0.6 and 1.3, respectively. 


3.2. Shock Heating 

In the shock heating scenario, hot electrons in the post-shock environment collision- 
ally heat pre-existing dust grains. Dwek (1987) and Dwek et al. (2008) provide a detailed 
description of this process for silicate dust grains and present post-shock equilibrium dust 
temperatures as a function of post-shock electron density, n e , and temperature, T e . Figure 
7 presents a similar analysis for graphite grains. 

For shock heating to occur, the forward shock must have sufficient time to reach the 
pre-existing dust grains. As discussed in further detail in Section 3.3, the peak supernova 
luminosity of SN 2005ip ( L peak ~ 10 9 L 0 ) vaporizes all dust grains out to a radius, r eva p ~ 
10 16 cm ( 0.01 ly). For the fastest observed shock velocities (v s ~15.000 km s~ x ), a shock 
would require ~70 days to reach the evaporation radius, which is consistent with the earliest 
observations of the HNI component (see Figure 6). The earliest observation of the WMI 
component does not occur until the Spitzer observations on day 936. by which point even 
the slower shocks would have crossed the evaporation radius. 

An independent calculation of the shocked dust mass provides a useful consistency check. 
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Assuming the hot, post-shocked gas heats the dust shell, estimates of the total gas mass will 
determine the dust-to-gas mass ratio. The upper limit on the volume of the emitting shell 
is given at 

When = 47rr s 2 Ar s , (5) 

where r s is the radius of the shock with velocity, v s , at an age, t, and A r s is the shell thickness 
defined by the distance traveled by the shock over the grain sputtering lifetime, r sputt , 

Ar shell = ^^s^sputt) (6) 

provided that r sputt < (t, r coo i), where r coo i is the radiative cooling time-scale. The factor of - 
comes from the shock jump conditions. For typical post-shock gas temperatures (> 10 6 K), 
Dwek & Arendt (1992) give the sputtering lifetime for a grain size, a, and gas density, n g , as 


Tsputt(yr) « 10 6 


a (/inn) 
n g ( cm -3 ) ’ 


The total gas mass of the emitting dust shell is therefore 


( 7 ) 


T/g ft-gWlllVs hell- 


( 8 ) 


Plugging in the above equations 

M g (M 0 ) « 8.3 x 10' 5 (: 


v s 


(IV, 

f a \ 

Vyr/ ' 

vgm/ 


. 1000 km s“ 

which reveals the mass of the shocked gas is independent of the grain density. 


( 9 ) 


Assuming a dust-to-gas mass ratio expected in the H-rich envelope of a massive star, 
z d=jf 0.01, gives the expected dust mass. For the maximum observed shock velocity, 
v s < 15,000 km s~\ and age, t = 936 days, Table 3 compares the predicted dust mass to 
the observed mass listed in Table 2. Assuming shock heating is also responsible for the 
HNI component at early times, Table 4 compares the predicted HNI dust mass at an age 
t = 70 days to the observed mass listed in Table 2. Although no measurement of the 
dust mass exists on day ~70, the relatively constant HNI luminosity (see Figure 6) and 
temperature (see Fox et al. (2009)) suggest a relatively constant HNI mass throughout the 
extent of the observations. 


For both the WMI component on day 936 and the HNI component on day 70, only 
large grains (a > 0.3 /mi) can reproduce the observed dust masses. These grain radii are 
large compared with typical grain sizes observed in supernova shocks (Dwek et al. 2008). 
Furthermore, upper limits were assumed for both the shock velocity (see Section 3.1) and 
dust-to-gas mass ratio (Williams et al. 2006). Lower values would require even larger grain 
sizes to reproduce the observed dust masses. These results likely rule out shock heating. 



3.3. Possible Emission From an IR Echo 


For an infrared echo scenario, the supernova luminosity heats a shell of dust at a radius, 
r, to a peak temperature, T d . This outer shell may be pre-existing at the time of the super- 
nova explosion or it may form when the peak supernova luminosity creates a vaporization 
cavity. In either case, light travel time effects cause the thermal radiation from the dust 
grains to reach the observer over an extended period of time, thereby forming an ‘IR echo 5 
(Bode k Evans 1980; Dwek 1983). The infrared luminosity plateau occurs on year long time 
scales, corresponding to the light travel time across the inner edge of the dust shell. As dust 
cools from the peak temperature, it will contribute flux at longer wavelengths. As noted in 
Section 2.3, however, the SN 2005ip spectrum is best fit by two components, as opposed to 
a continuous temperature distribution. Therefore, this analysis assumes a simple light echo 
model that is dominated by flux from only the warmest dust with a single temperature, T d . 

The equilibrium dust temperature is set by balancing the energy absorbed and emitted 
by the dust grains, 

Tabs = T ra( j, (1H) 


where, for a single dust grain, 


and 


r 3 2 7102 

L abs - 47rr SN 47rr 2 

J 7rR l/ (T SN )Qabs(i / )^ 

^bol 7T& i 

<TT s 4 n 4 r*J 

^ Ri'(TsN)f?abs(l / )^^ 

Trad = 47TU 2 J 

7T B v (Td)Qabs(^ )du 

16 o 

= T «pa 

J B v {T d )K(u)dv 


( 11 ) 


( 12 ) 


where tsn is the effective supernova emitting radius, T bo i is the UV-optical luminosity, and 
Tsn is the effective supernova blackbody temperature. T bo i follows from Equations 10, 11, 
and 12, 


T b oi — 


64 2 4 f B u (T d )K(v)dv 

3 Pa ' a SK f 


(13) 


Tboi depends on the grain radius because although the dust opacity coefficient, k, in Figure 
4 is independent of grain radius for thermal emission at longer wavelengths, it does depend 
on grain radius for absorption at shorter wavelengths (e.g., UV and optical). 


Using equation 13, Figure 8 plots contours of the shell size, r, as a function of both 
luminosity, L bo i, and observed dust temperature. T d . The luminosity is treated as a central 



point source, assuming the emitting region is internal to a spherically symmetric dust shell. 
Several grain sizes are considered for dust with a graphite composition. Although the cal- 
culation assumes Tsn ~ 10,000 K, the result is fairly insensitive to this choice. The vertical 
lines show the observed graphite dust temperatures listed in Table 2 and the approximate 
vaporization temperature of graphite dust, T ev ap ~ 2000 K. The horizontal lines show both 
the observed peak, L pea k» and late- time optical/infrared plateau, Lpiateau, luminosities from 
Figure 6. 

The shaded regions highlight the dust shell radii allowed by the constraints. The shock 
radius sets the lower limit as any dust within this radius would be independently heated 
or destroyed by the forward shock (see Section 3.2 above). (For the WMI component, the 
blackbody radius, r hh (W MI), actually sets the minimum radius, as described in Section 3.1.) 
Both the fast, r sl , and slow, r s2 , shocks described in Section 3.1 are considered, distinguished 
by the hashed region. The near-infrared plateau time-scale sets the upper limit. Figure 6 
shows the plateau extends for at least ~2.6 years and still shows little sign of declining. The 
current upper limit is therefore set by a “plateau radius,” r p , of 1.3 light years, although 
the final upper limit will only be determined once the light-curve plateau begins to decline. 
Should the plateau extend even longer than 2.6 years the larger implied radius only makes 
a light echo less likely. 

Three possible light echo scenarios exist: 

1) In the first scenario, a sphere of dust completely encompasses the progenitor at the time of 
the explosion. The peak luminosity vaporizes all dust within a radius, r eV ap, and warms the 
inside of the remaining dust shell to nearly the vaporization temperature (T evap ~ 2000 K). 
Figure 8(a) shows that the observed peak optical luminosity, L pea k, yields a vaporization 
radius of r evap ~ 0.013 light years for a = 0.01 /.an graphite grains and an even smaller 
radius for larger grains (see Figures 8(b) and 8(c)). These small evaporation radii, however, 
are inconsistent with the observations as a shell of this size cannot produce an infrared echo 
on 3 year time scales. Furthermore, these evaporation radii are smaller than both shock radii. 
Even if the actual peak luminosity were a factor of 5 larger than observed (L est ~ 5 x 10 9 L e ), 
the vaporization radius would be insufficient to account for the observations. We therefore 
rule out this scenario for both the HNI and WMI components. 

2) The dust shell inner limit need not lie exactly at the vaporization radius. If the 
progenitor underwent an eruption many years before the supernova, the dust shell may lie 
at larger radii. In this second scenario, the peak luminosity, L pea k, heats the dust shell 
inner radius to only the observed temperature. The light echo duration therefore defines 
the minimum cavity radius (i.e., r p — 1.3 light years). Figure 8 shows that a minimum peak 
luminosity of Lboi > 5 x 10 10 Lq is required to heat a dust shell of 1.3 light years in radius to 



the observed WMI temperature, while a minimum peak luminosity of Lboi > 5 x 10 11 L 0 is 
required for the HNI component. A larger shell radius requires even larger peak luminosities. 

These required peak luminosities are significantly larger than the observed peak lumi- 
nosity. The true peak luminosity of SN 2005ip, however, is not well constrained. The earliest 
R-band photometry was obtained 14-days post-discovery, and the discovery occurred a few 
weeks following the actual explosion so that the peak luminosity was likely several times 
larger than the observed early-time Abband luminosity (L e st ~ 5 x 10 9 L 0 ). Still, while a 
significant amount of optical absorption might be expected by large amounts of pre-existing 
dust, significant reddening was not observed (Smith et al. 2009b) and a peak luminosity 
> 10 10 L 0 is unlikely from this extrapolation. 

The shock breakout in the minutes to hours following the supernova explosion may reach 
peak luminosities > 10 11 L 0 (Soderberg et al. 2008; Rest et al. 2009), but no such breakout 
was observed for SN 2005ip. A peak luminosity > 10 11 L 0 would have made SN 2005ip one of 
the most luminous core-collapse events ever observed (Quimby et al. 2007; Rest et al. 2009), 
but figure 11 of Rest et al. (2009) suggests SN 2005ip is nearly an order of magnitude fainter 
than SNe 2006gy 2008es, 2005ap, and 2003ma, especially at early times. Furthermore, 
optical emission from the late-time circumstellar interaction successfully accounts for the 
observed dust temperatures (see below). These reasons rule out a light echo driven by the 
peak supernova luminosity for both the HNI and WMI components. 

3) A final scenario considers a pre-existing dust shell similar to scenario 2 above, but in 
this case the shell’s inner radius is located at an intermediate radius between the shock and 
plateau radii (as defined on day 936). The late-time optical emission, I p i a teau, continuously 
heats the dust shell to the observed temperature. This scenario is not so much a traditional 
infrared echo as it is a reprocessing of the optical emission by the dust. (Some authors refer 
to this as a circumstellar shock echo (Gerardy et al. 2002).) If the circumstellar interaction 
occurs on a time scale greater than the light travel time across the dust shell, the shell radius 
does not set the infrared plateau length. The observed flux therefore accounts for the entire 
shell. 

Figures 8(a) and 8(b) show that for both a = 0.01 and 0.1 pm graphite grains, L p !ateau 
can heat a dust shell of radius r « 0.01 ly to Thni and a shell of radius r % 0.05 ly to Twmi- 
For a = 0.5 /an graphite grains, Figure 8(c) shows that L p i ateau can only heat a dust shell of 
radius r < 0.005 ly to T H ni and a shell of radius r < 0.04 ly to Twmi. fa this scenario, both 
the fast {r S \ = 0.125 ly) and slow (rsa = 0.042 ly) shock radii are larger than the HNI shell 
radii by day 936, ruling out a light echo of this sort for the HNI component. The same is 
true for the WMI shell composed of larger grains (a > 0.5 pm). In the case of the WMI shell 
composed of smaller grains (a < 0.1 pm), however, the slower shock has not yet reached the 
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shell radius (r « 0.05 ly). Not only is this scenario possible for the WMI component, but the 
radius is consistent with the WMI blackbody radius (r bb (WMI) = 0.048 ly). Furthermore, 
the grain sizes are typical of those observed in other supernova circumstellar environments 
(Dwek et al. 2008). 

Predictions for the light curve evolution can be made from this model. Dust that 
remains at radii beyond the slower shock radius, r s 2 , will continue to radiate and contribute 
to the light echo plateau. Although the exact distribution of these dense, dusty regions 
is not well-known, this scenario suggests they must be distributed in such a way that the 
fastest shocks do not interact with the dust. As described in Section 3.1, these dense regions 
may be concentrated in either an equatorial disk or clumps. The clumps must have a large 
filling factor if the fastest shocks are not to interact with a significant portion of the dust. 
As the slower shocks continues to expand, however, they will ultimately destroy the dust. 
Assuming the dust lies at a radius consistent with 0.1 pm grains (see parameters in Table 2), 
the WMI flux will begin to decrease at t « PALilil ~ 2775 days post-discovery and continue 
to decrease as a function of the emitting shock radius and the dust distribution. While this 
scenario is entirely consistent with the WMI observations at the current time, mid-infrared 
observations at later epochs can reveal the accuracy of this model’s predictions. 


4. Discussion 
4.1. The HNI Component 

The above HNI component analysis rules out the possibility of both pre-existing dust 
scenarios (i.e., shock heating or infrared echo). Only condensation models (i.e., in the ejecta 
or cool, dense shell) remain as viable scenarios for explaining the origin of the HNI compo- 
nent. Smith et al. (2009b) confirm new dust formation via extinction in the broad (~15,000 
km s -1 ) Ha wings between days ~60-170 and in the intermediate (~2000 km s -1 ) He I wings 
at days >413, but the location of this new dust remains ambiguous. The extinction in the 
broad and intermediate components suggest new dust in the fast ejecta at early times and 
post-shock cool, dense shell at later times, respectively. 

To form dust at the high ejecta velocities, however, is difficult due to the lack of heavy 
metals traveling at these speeds. If the dust formed in the cool, dense shell at early times, as 
suggested in the case of SN 2006jc. (Mattila et al. 2008), this would explain the attenuation of 
the broad Ha line observed by Smith et al. (2009b) without having to invoke dust formation 
at high velocities, but would fail to explain the lack of observed intermediate width He I 
lines produced by the circumstellar interaction. Smith et al. (2009b) propose one alternative 
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scenario in which the dust may form in the post-shock gas of individual clumps, which then 
become incorporated into the expanding fast ejecta when a clump is eventually destroyed. 
This scenario, however, requires several assumptions, most of which rely on unknown clump- 
ing properties. While the available observations limit the ability to isolate the precise region 
of dust formation, the relatively flat near-infrared flux in Figure 6 suggests a majority of the 
dust contributing to the HNI flux must have formed at early epochs. 

Newly formed dust cannot reproduce the observed HNI luminosity plateau solely by 
cooling from condensation to the observed temperature T H ni ~ 800 K. Given an average 
energy per particle e = C g ATd for the specific heat for graphite C g (given by Draine & Li 
(2001)) and AT ~ 1200 K, an unsustainable mass condensation rate of M = ThniA ~ 100 
M 0 day -1 would be necessary to reproduce the observed flux (see a more detailed explanation 
of this argument for the case of SN 2006jc in Fox et al. (2009)). 

Instead, an alternative heating mechanism must power the thermal emission from this 
newly formed dust. Figure 6 shows that radioactive heating is insufficient to power the 
late-time near-infrared emission. While heating by the reverse shock is possible, a similar 
analysis as performed in Section 3.2 suggests this scenario is unlikely because only large 
grains (a > 0.5 /.tm) can reproduce the observed flux. More likely, the optical luminosity 
generated from the forward shock interaction continuously heats the newly condensed dust 
in the same way that it heats the WMI component discussed in Section 3.3. Unlike the 
infrared echo scenario, however, the newly formed dust exists interior to the forward shock 
radius and the shock emission cannot be treated as a central point source. As a consistency 
check, the comparable near-infrared and optical luminosities in Figure 6 yield an optical 
depth, r ss £ n1 r+l 0 pt = 0-5, which is comparable to the the optical depth given for the HNI 
component by equation 4 in Section 3.1. 


4.2. A Multi-Component Model 

A multi-component dust model for SN 2005ip now begins to emerge, composed of both 
an inner, “hot” (~800 K), near-infrared (HNI) and an outer, “warm” (~400 K), mid- 
infrared (WMI) component. Newly formed dust in the either the ejecta or cool, dense 
shell likely dominates the HNI component. The WMI temperature and blaekbody radius, 
rjb(WMI) «4.8 x 10 1 * cm (0.048 ly), are consistent with a pre-existing dust shell heated by 
the observed late-time optical luminosity generated by the forward shock interaction (Sec- 
tion 3.3). Figure 9 illustrates the locations of each component, as well as the likely origins 
and heating mechanisms. 
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The large pre-existing dust mass that contributes to the WMI flux component suggests 
significant mass loss from the progenitor. Assuming a dust-to-gas ratio Z d = 0.01. the 
observed WMI dust mass listed in Table 2 (M<f(WMI) ~ 0.01 — 0.05 M 0 ) yields a total gas 
mass of M g (WMI) ~ 1 - 5 M 0 . This mass accounts for the entire WMI emission, as the 
entire shell contributes to the observed flux given the size of the shell (r^fWMI) ~4.8 x 
10 16 cm (0.048 ly)) is significantly less than the observational time-scale (936 days). The 
associated mass loss rate is 


M 


M g (WMI) 


A r 

7.5 x 10 


-3 ( Mg (WMI) \ / 

/V 120 km 


M, 


© 


:)(^)( 


A rJ 


dvi© 


yr- 1 ). 


(14) 

(15) 


On day 413 post-discovery, Smith et al. (2009b) measured a progenitor wind velocity v w = 
120 km s -1 . Assuming a constant wind velocity and thin shell = A), the mass loss 
rate is M « 7.5 x 10“ 2 — 3.8 x 10 _1 M 0 yr -1 , which is two-three orders of magnitude 
larger than calculated by Smith et al. (2009b) at inner radii. Furthermore, this mass loss 
rate is likely a lower limit. He I P-Cygni profiles at 1.083 /im from the TripleSpec spectra 
show the progenitor wind velocity may have been closer to 200 km s -1 at the WMI radius 
(Figure 10), assuming the P-Cygni feature is generated by the cool, low velocity circumstellar 
environment that coincides with the WMI shell. 


The larger mass loss rate (and possibly faster progenitor wind velocity) of the WMI shell 
suggests a denser progenitor wind (or “eruption”) occurred at t < -■■■-■A.Mi) ~ ~ 

125 years prior to the core collapse. Although Smith et al. (2009b) conclude that the SN 
2005ip progenitor was likely a Red Super Giant (RSG), these stars typically only have wind 
speeds v w ~ 20 — 40 km s -1 and mass loss rates up to M = 10~ 4 — 10~ 3 M© yr -1 (Smith et al. 
2009a). The observed characteristics associated with the WMI shell are more consistent with 
Luminous Blue Variable (LBV) stars (e.g., Davidson 1989; Humphreys & Davidson 1994), 
which can have wind speeds on order of hundreds km s -1 (e.g., Leitherer 1997; Kotak k 
Vink 2006) and can have mass loss rates up to M = 10 _1 M© yr -1 (Smith k Owocki 2006; 
Smith k Hartigan 2006; Smith et al. 2007). 


5. Conclusion 

The Spitzer spectrum presented here confirms the presence of warm dust in SN 2005ip. 
Combined with near-infrared observations, the results show evidence for two independent 
dust masses: a hot, near-infrared (HNI) and warm, mid-infrared (WMI) component. Infrared 
observations span the peak of the thermal emission, thereby providing strong constraints on 
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the dust, mass, temperature, and luminosity, which serve as critical diagnostics for disentan- 
gling the origin and heating mechanism of each component. The HNI dust mass originates 
primarily from newly formed dust in the ejecta, while the WMI component likely origi- 
nates from an circumstellar shock echo that forms from the heating of a large, pre-existing 
dust shell. In both cases, the heating mechanism is likely the optical luminosity gener- 
ated from the forward shock interaction with the circumstellar medium. For wind speeds of 
~ 100 — 200 km s“\ the WMI dust shell likely formed via an eruption ~ 100 years before 
the supernova. These characteristics are consistent with a LBV progenitor, which has been 
linked to some core-collapse supernovae (Kotak & Vink 2006), but is an emerging trend par- 
ticularly within the Type Iln subclass (e.g., Gal- Yam et al. 2007; Smith et al. 2007; Trundle 
et al. 2008; Smith et al. 2008a, 2010; Gal- Yam & Leonard 2009). 

The newly formed HNI dust mass, M d (HNI) ~ 5 x 10”~ 4 M 0 , is at least two orders 
of magnitude lower than predicted by ejecta condensation models to reproduce the large 
amounts of dust observed at high redshifts (Todini & Ferrara 2001; Nozawa et al. 2003, 
2008). A dust mass of ~ 10~ 4 M e is comparable to other Type Iln events (e.g., Gerardy 
et al. 2002; Pozzo et al. 2004; Meikle et al. 2007). For clumpy ejecta, it should be noted 
that both semi-analytical models (Varosi & Dwek 1999) and Monte Carlo radiative-transfer 
simulations (Ercolano et al. 2007) have shown that dust masses can be at least an order of 
magnitude larger than that predicted by the analysis presented in Section 2.3. The clump 
structure for any supernova, however, remains unconstrained at present. The pre-existing 
WMI dust mass, Afd(WMI) ~ 0.05 M 0 , is quite significant. Although, in the case of SN 
2005ip, the forward shock may ultimately destroy this dust (mid-infrared observations at 
later epochs will reveal the accuracy of this model), the forward shocks may be sufficiently 
decelerated in other systems to allow for dust survival. 

Late-time mid-infrared observations of dust in the supernova environment provide for a 
unique interpretation of the circumstellar environment and progenitor system. Nonetheless, 
mid-infrared observations are quite rare. Presented in this paper is the first mid-infrared 
spectrum of any Type Iln supernova. In the future, we hope to grow the database of mid- 
infrared observations of Type Iln supernovae. Doing so, however, can be slow as Type Iln 
events are particularly rare, consisting of only ~2-3% of all core-collapse supernovae (Gal- 
Yam et al. 2007) and occurring at a rate of no more than 10/yr out to 150 Mpc (Dahlen & 
Fransson 1999). We therefore plan to revisit all Type Iln supernovae from the past ten years 
to determine the degree to which this subclass exhibits late-time dust emission and identify 
the emission mechanism. 


This work is based on observations made with the Spitzer Space Telescope (PID 50256) 
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Table 1. Spitzer Observations Summary 


Target 

Field 

Position 

RA 

DEC 

Flux 

Density 

Bands 

Int/ 

Pixel 

(secs) 

Total 

Duration 

(secs) 

Ramp x 
Cycles 
(sec x #) 

# of 
AORS 

SN 2005ip 

09:32:06 

+08:26:44 

0.4 mJy 
(8/im) 

IRAC 

all 

300 

1062 

- 

1 

SN 2005ip 

09:32:06 

+08:26:44 

0.4 mJy 
(8pm) 

IRS 

Spectral Mapping 
low short both 

- 

9516 

60 x 12 

1 
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Table 2. Fitting Parameters 


Fit 


1.0 /nr i(warm) + 0.3 /an (hot) 

0.3 /xm(warm) + 0.3 /mi (hot) 

0.1 /tm(warm) + 0.3 /an (hot) 
0.01 /nn(warm) + 0.3 n m(hot) 
0.001 /jm(warm) + 0.3 /an (hot) 

1.0 /tm(warm) +0.1 /an (hot) 

0.5 /im(warm) +0.1 /xm(hot) 

0.3 /an (warm) +0.1 /an (hot) 

0.1 /an (warm) +0.1 /rm(hot) 
0.01 /nn(warm) +0.1 /mr(hot) 
0.001 /im(warm) + 0.1 /an (hot) 

1.0 /an (warm) + 0.01 /an (hot) 
0.3 /xm(warm) + 0.01 /an (hot) 
0.1 /an (warm) + 0.01 /an (hot) 
0.01 /tm(warm) + 0.01 /tm(hot) 
0.001 /a n(warm) + 0.01 /tm(hot) 
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Table 3. Shock Heating Mass Predictions at t — 936 Days for v s = 15,000 km s 1 



WMI 

HNI 

a (/an) 

M d (M 0 ) 
Predicted 

M d (M©) 
Observed 

M d (M©) 
Predicted 

M d (M 0 ) 
Observed 

0.01 

1.8e-4 

4.5e-2 

1.8e-4 

5.7e-4 

0.1 

1.8e-3 

4.3e-2 

1.8e-3 

5.2e-4 

0.3 

5.5e-3 

2.8e-2 

5.5e-3 

4.1e-5 

0.5 

9.2e-3 

1.2e-2 

- 

- 

1.0 

1.8e-2 

5.9e-3 

- 

- 
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Table 4. Shock Heating Mass Predictions at t = 70 Days for v s = 15,000 km s 1 



HNI 

a (/Jin) 

M d (M 0 ) 

M d (Mo) 


Predicted 

Observed 

0.01 

1.0e-6 

5.7e-4 

0.1 

1.0e-5 

5.2e-4 

0.3 

3.1e-5 

4.1e-5 



Fig. 1. — Spitzerj IRAC and IRS observations of SN 2005ip in NGC 2906 obtained in June, 
2008 (^-943 and 936 days post-discovery, respectively). Shown here is the Post-BCD 3.5 
//m IRAC image with the IRS map overlay. SN 2005ip is directly above the galaxy, while 
the galaxy nucleus also falls within the mapping scheme. 
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Pixel Number 


Fig. 2. — Spatial profile showing SN 2005ip centered in the slit with the galaxy nucleus at 
the left shown by the AdOpt tool within the SMART data analysis package. The two point 
sources are fit simultaneously along with the underlying background. 




Fig. 3.— Combined APO/TripIeSpec near-infrared spectra and Spitzer/IRAC and IRS mid- 
infrared data. Although the epochs do not match perfectly, TripleSpec observations on day 
1243 post-discovery (not shown) suggest little evolution occurred since the day 862 spectrum 
plotted here. We therefore treat the combined spectra as a single spectrum. The thermal 
emission that dominates the spectrum confirms the presence of warm dust. 



k (Graphite) 
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Q (Graphite) 


Q (Silicate) 


Fig. 4. — The dust opacity coefficient, k, and emission efficiency, Q, for both graphite and 
silicate for several grain sizes as a function of wavelength. The emission efficiency is given by 
the equation 3. The dust opacity coefficient is nearly independent of grain radius assuming 
grains with radii less than the emitting wavelengths (i.e. a < 1 pm at infrared wavelengths). 
For large grains (a > 1 /an), the dust can be approximated by a blackbody (Q = 1). 
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X (pm) 

Fig. 5. — Best fits of the graphite and silicate models given by equation 2 to the combined 
near- and mid-infrared data. The lack of an emission feature at ~9 /mi immediately rules out 
the silicate model. The combined spectrum is best fit by a multi-component model, consisting 
of both a "hot” (~800 K), near-infrared (HNI) and “warm” (~400 K), mid-infrared (WMI) 
component. 
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Fig. 6. — The evolution of SN 2005ip since discovery. Included in this plot are near-infrared 
photometry from Foxet al. (2009), unfiltered (R- band) photometry from Smith et al. (2009b), 
and the WMI and HNI luminosities on day 936 post-discovery derived from the best fits in 
Section 2.3 and shown in Figure 5. The radioactive decay for a typical (0.1 M© of o6 Ni) 
supernova is overplotted (dashed line). 



T, 

Fig. 7. — Post-shock equilibrium dust temper? 
electron density, n e , and temperature, T e . Fo: 
T e % 10 : K. the grains become increasingly tra 
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Fig. 8.— Contours of the dust shell size, r, plotted as a function of both the supernova lumi- 
nosity, Lboh and dust temperature, Td- given by equation 13. Several grain sizes are consid- 
ered for dust with a graphite composition. The vertical lines show the observed graphite dust 
temperatures listed in Table 2 and the approximate vaporization temperature of graphite 
dust. r evap ~ 2000 K. The horizontal lines show both the observed peak. T peak , and late-time 
optical plateau, I P s a teau- luminosities from Figure 6. The shaded regions highlight the dust 
shell radii allowed by the constraints. The shock radius sets the lower limit as any dust 
within this radius would be independently heated or destroyed by the forward shock f,see 
Section 3.2). Both shock radii described in Section 3.1 are considered. The horizontal line fill 
distinguishes r s L from r s 2 . The near-infrared plateau time-scale sets the upper limit. Figure 
6 shows the plateau extends for at least ~2.6 years and still shows little sign of declining. 
The minimum upper limit is therefore 1.3 light years, although it may certainly be larger. 
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Fig. 9. — Illustration of the proposed multi-component model for SN 2005ip consisting of 
both a hot, inner (HNI) component and warm, outer (WMI) component. The HNI dust 
mass originates primarily from newly formed dust in the ejecta, while the WMI component 
likely originates from an circumstellar shock echo that forms from the heating of a large, pre- 
existing dust shell. Collisional radiation from continuous shock interaction with the dense 
circumstellar medium generates the optical luminosity plateau (Smith et al. 2009b), which 
is the likely heating mechanism for both the HNI and WMI components. 
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Fig. 10. — The 1.083 /mi He I line. The shocked gas creates the broad (~ 1000 km s -1 ) emis- 
sion feature. The unshocked, slow moving, circumstellar medium formed by the progenitor 
wind produces the superimposed P Cygni profile. The blue edge of the P Cygni absorption 
feature. ~200 km s -1 , measures the wind velocity on clay 862 post- discovery. 





